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ABSTRACT 


Orthogonality  relations  between  the  eigenvectors  of  damped  linear 
dynamic  systems  with  lumped  parameters  are  derived;  and  from  these  relations 
coordinates  are  found  in  terms  of  which  uncoupled  equations  of  motion  can  be 
written.  Methods  are  developed  for  determining  transient  stresses  in  terms  of 
these  coordinates.  The  present  treatment  is  extended  to  systems  involving 
transient  damping  and  to  continuous  systems. 
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INTRODUCTION 


The  equations  of  motion  of  a  damped  linear  dynamic  system  with 
lumped  parameters  can  be  written  in  matrix  notation  as 

[%mj  hut,)  +va]fft}  =  mi 

The  classical  method  of  solving  these  equations  is  to  find  the  normal  modes  of 
oscillation  of  the  homogeneous  equation 


HM1  +MM  =  til 

and  from  these  modes  to  determine  a  set  of  normal  coordinates  in  terms  of 
which  there  is  no  inertial  or  elastic  coupling  between  the  equations  of  motion. 
However,  unless  the  matrix  jr]  happens  to  be  proportional  to  either  [m]  or 
[k],  velocity  coupling  still  exists,  and  the  chief  object  of  using  normal  coord¬ 
inates  is  defeated  except  when  the  coupling  terms  are  small  and  can  be  neglected. 

The  purpose  of  the  present  paper  is  to  derive  orthogonality  relations 
between  the  homogeneous  solutions  of  Eq.  (1),  and  from  these  relations  to 
develop  a  set  of  coordinates  in  terms  of  which  completely  uncoupled  equations 


of  motion  can  be  written. 


HOMOGENEOUS  SOLUTIONS 


Equation  (1)  can  be  written  in  reduced 


1 


form  as 


[R]{Z]  +  [/C]  (z  J  =  {Fit)] 


where 


[R]  = 


loj 

Ml 

M  ~ 

W  J 

Co] 


j 


To  obtain  the  homogeneous  solution  of  Eq.  (3),  let 


(z(t)  ]  = 


This  gives 


«[*]{*]  +  [*]{$]  =  { 0 ] 


(3) 


(4) 


(5) 


Equation  (5)  may  be  manipulated  into  a  form  more  amenable  to  matrix  iteration, 
such  as 


[um  =  i  w 


"Elementary  Matrices",  R.  A.  Frazer,  V/.  J.  Duncan  and  A.  R.  Collar, 
Cambridge  University  Press,  London,  1946,  p.  327. 
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where 


[Uj  =  -  [#]~'[R J 


r  m  ru 

-[cJM  ~[cl[r] 

L  J 


and  where 

W  =  [Jr] 


The  matrix  [c]  is  a  set  of  flexibility  influence  coefficients. 

The  solution  of  Eqs.  (5)  or  (6)  will  yield  2N  eigenvalues  and 
eigenvectors. 


anti 

n 


ns  !,  Z,  •  ■  •,  IN 


where  N  is  the  number  of  degrees  of  freedom.  For  a  stable  system  each  an  is 
either  real  and  negative  or  complex  with  a  negative  real  part.  The  complex 
roots  occur  as  complex  conjugate  pairs  with  corresponding  complex  conjugate 
modal  columns.  A  relatively  simple  method  of  iterating  Eq.  (6)  for  complex 
eigenvalues  and  eigenvectors  is  described  by  Frazer,  Duncan  and  Collar  2. 

If  the  matrices  [m],  [r]  and  [k]  are  symmetric,  the  matrices  [R]  and 
[K]  are  symmetric,  and  Eq.  (5)  can  be  written  in  either  of  the  two  ways 

2  Ibid 
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«»[*]{§"]  +  IK]  ($(n))  ={ 0 } 


(7) 


and 

j  +  = (°v 

where  the  superscript  T  denotes  a  transposed  matrix.  The  premultiplication  of 
Eq.  (7)  by  ($^J  T  and  the  postmultiplication  of  Eq.  (8)  by  yields 

+  {i'«}T0<]{iw}  =  ° 

and 


Jf#wJ  +{iMfd<mn  =  0 

When  Eq.  (10)  is  subtracted  from  Eq.  (9),  one  obtains 

K-  «„)($MiT[*3teM}  =  o 

Thus,  unless  a  =a  ,  the  orthogonality  relations  are 
'  m  n  a  / 

{$lm^]T[f?]{$tn)]  =  0  ,  when  m  *  n 

and 

~  0  ,  when  m  &  n 


(10) 


01) 


02) 


(13) 
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where 


R„  =(iM}Tm{$“"}  =  +  i^Vtr W} 


and 

=  {i*>}TlFtt)]  = 


From  the  orthogonality  relations  it  can  also  be  shown  that 

(m}  =  10 

h  *r\ 

and 


l/m)  =  |2!  [«.  M  +  [rj]  {*<■”>]  (19) 

also  that 

=  (20) 

and 


I  ^  f/*J  =  foj 

n  Kn 


(21) 


For  zero  initial  conditions,  the  solution  of  Eq.  (17)  is 


(22) 
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DYNAMIC  STRESSES 


^  "Methods  in  Transient  Stress  Analysis",  R.  L.  Bisplinghoff,  G.  Isakson  and 
T.  H.  H.  Rian,  Journal  of  the  Aeronautical  Sciences,  Vol.  17,  No.  5, 
p.  259,  May  1950. 
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where  the  double  star  indicates  that  only  "rigid  body"  modes  are  included  in 
the  summation.  Therefore, 


rn(t) 

-  <Xy>  Rn 


(33) 


and  so. 


O'. 

‘  (static) 


=  Ls'/)j  {/(*)]  -  Z  K 


n 


(34) 


Thus,  through  the  use  of  Eq.  (31),  the  stress  due  to  the  applied  forces  can  be 
accounted  for  exactly  even  when  a  limited  number  of  vibration  modes  are  used. 


TRANSIENT  DAMPING 

Many  problems  in  aeroelasticity  require  that  the  velocity  in  Eq.  (1) 
be  replaced  by  a  convolution  integral  containing  the  Wagner  function;  i.e.,  a 
single-degree-of-freedom  equation  of  the  same  type  as  Eq.  (1)  would  become 


m 


t 

f  +  rj  <f(r)  +  k  <f  =  fit) 


(35) 


where  4>(t)  is  the  Wagner  function  which  accounts  for  unsteady  flow  effects. 
It  can  be  shown  by  means  of  the  Laplace  transform  that,  if  the  usual  exponential 
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approximation  of  the  Wagner  function  is  used,  Eq.  (35)  is  equivalent  to  a  fourth 
order  differential  equation  of  the  type 

aT  +  *>?  -hc$  +  di+e<i^  g(t) 

Thus,  Eq.  (1)  becomes 


(36) 


[%]{%}  +f +fc;][%1  +  fcjftj}  =  ho] 


(37) 


This  can  also  be  expressed  in  reduced  form  as 


(38) 


where  now 


[0]  [0]  [of 

~[oj  [o]  -w  roj~ 

[o]  [o]  M  W 

M  = 

[oj  -w  -w  roj 

[0]  M  W  M 

-w  -JM  -fcJ  fetf 

M  W  CO  W 

roj  r^j  i>j 

1 

M 

r-^-N 

f*J 

,  (Fm]  = 

ft/ 

M. 

ASRL  TR  25-20 


11 


CONTINUOUS  SYSTEMS 


Orthogonality  relations  for  continuous  systems  can  be  obtained  in  a 
manner  analogous  to  that  for  lumped  parameter  systems.  Consider  for  example 
the  damped  motion  of  a  simple  beam.  The  equation  of  motion  is 


[*' 

i  (u>  t)  +  r(y)  z  (y,  t)+  k(H,  l)  z  (y „ 1)  M  y  z=  /  ft/  t) 

where  k(y,v)  is  a  stiffness  influence  function^.  This  equation  can  be  written 
in  reduced  form  as 


(39) 


m(!/)±(y,t)  -  m(j)  K  (*,t)  =  0 

»(y)  ffy/t)  +  rd)  z[i,t)  +  Jk(!,'j)  z(9,t)  <*1  -  J 


(40) 


where 


9. 


To  obtain  the  homogeneous  solution  of  Eqs.  (39)  or  (40),  let 


ziy.t)  -  4>(y)  e**  an<4  ? 6o *)  =  &(*)  e 


or  t 


"Aeroelasticity",  R.  L.  Bisplinghoff,  H.  Ashley  and  R.  L.  Halfman, 
Addison-Wesley  Publishing  Company,  Inc.,  Cambridge,  Massachusetts, 

p.  26. 
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Equation  (40)  becomes 
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In  Eq.  (43)  the  subscripts  m  and  n  may  be  interchanged  to  give 


°<m  J (*(/)  <£„[*) /3„(y)  +  r(y)  fay)  fay)  +  m(y)  /3„(y)  fay))  dy 


+ j  ^  j k  ^m(v)  <P„  id)  d  9  -  m(y)  /in[y) j  dy  -  0 


When  Eq.  (44)  is  subtracted  from  Eq.  (43),  one  obtains 


(44) 


K(»)  /**,(*)  +  r(y)  tj?)  +  m0')'3Jy)  t„(y)) 


o 


(45) 


because  k(y,y)  is  a  symmetric  kernel.  Thus,  the  orthogonality  relations 
are 


/ (mto)  tnto) /**,(#)  +  ’iy)  4>M  fa(y)  +  m(y)/3„(y)  ^[y)J dy  -  o/  wt,*,, 


nt  yz  n 


(46) 


and 


k  (*'  i)  ^„to)  4>m(d)  dy  -  m(d)  /3»(y)  /$„(?))  d]  =  0,  when  w>  dn 


(47) 


Through  these  results  uncoupled  equations  of  motion  similar  to 
Eq.  (17)  can  be  derived  for  this  continuous  system. 
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